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APPENDIX I 



% Creating a network topology obj 

network__topo = topo('init'); 

addlink(network_topo); 

labelnames(network_topo); 

save network__topo; 



% graphically place nodes on screen 
% graphically connect up nodes 
% graphically label nodes 

% save network_topo for future use 



% Top level procedure to compute paths that optimize use of network capacity 
% inputs: 

% D = traffic demand matrix 

% (retrieved from predictions stored in TMS Statistics Repository) 
% network_topo = topology object defining the network topology 
% P = network policy information 

% (matrix of reserved capacity, which indicates links whose use 

% is administratively prohibited or which should not be 

% completely allocated) 

% outputs: 

% allocated_paths() = list of paths to set up, to TMS signalling system 



C = capacity(network_topo); % retrieve network topology information 

C = C-P; 

saved_C - []; 
saved_SLA = []; 
assigned jaths = []; 
round = 0; 

[SLA, S] ^ create_ordered_sla(D); 

F-SLA(l) 

for F- SLA, 

rounds round +1; 
saved_C {round} = C; 
saved_SLA{round} =^ F; 

F % display the flow 

W = calc_weights('calcweight2',F,C); 
[dist, P] = floyd(W); 
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path = findpath(P,F.i,F j); 

assigned_paths {round}. path = path; 
assigned_paths {round}. flow ^ F; 

if (isempty(path)) 

fprintf(i;no path for flow:\n'); F 

else 

C = compute_residual_capacity('c - F.bw',path,F,C); 

end 

end 

function [W] = calc_weights(func,F,C) 
% function [W] calc_weights(func,F,C) 
% 

% Compute the weights by calhng func on each elt of C 

% func must be of the form double func(Flow F, Capacity_elt c, node i, node j) 

func = fcnchk(func); 

for i = l:size(CJ) 

for j = l:size(C,2) 

W(ij)-feval(func,F,C(ij)4j); 

end 

end 

function [w] = calcweight2(F,c,ij) 
% function [w] = calcweight2(F,c4 j) 
% basic weight calc 

if(0 = c) 

w = inf; 
return; 

end 



% rule out paths that can't hack it 
if(F.bw>c) 

w = inf; 

return; 

end 
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w = 1 / (c - F.bw) ; % fill links with most capacity first 

function [C] = compute_residual__capacity(fixnc, path, F, C) 
% function [C] = compute_residuaLcapacity(func, path, F, C) 

% 

% Update capacity characteristics in C to reflect flow F being 

% allocated along path using function func 

% func should be of the form 

% C_element func(C_element c, Flow F) 



if(length(path) 1) 
return; 

end 

func =^ fcnchk(func;cVF'); 



index = 1 ; 
src = path(index); 
20 index = index + 1; 

for index = index :length(path) 
dst path(index); 

25 C(src,dst) = feval(func,C(src,dst),F); 

src = dst; 

end 

30 function [SLA, S] = create_ordered_sla(D) 

% function [SLA] = create_ordered_sla(D) 
% takes the demand matrix and returns a list of SLAs, 
% SLA of the form [ struct ; struct ; ... ] where struct is [BW, i, j] 
% S of the form [ [BW, i, j] ; [BW, i, j] ; ...] 



s-[]; 



for i = l:size(D,l) 
40 forj = l:size(D,2) 

if(D(ij)-= 0) 

S = [[D(ij)ij];S]; 

end 

end 

45 end 
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[Y,I] = sortrows(S,l); 

S = Y(size(Y,l):-l :!,:); % reverse order 

5 

SLA = structCbw',num2cell(S(:,l));i',num2cell(S(:,2));j',num2cell(S(:,3))); 
return; 

10 function [path] = findpath(P,i,j) 

% function [path] = findpath(P) 
% 
% 

15 path=[]; 

if(i==j) 

path= [i]; 
return; 

20 end 

if(0 = P(i,j)) 

path = []; 

else 

25 path = [findpath(P,i,P(iJ)) j]; 

end 

function [D, P] = floyd(W) 
% function [D, P] = floyd(W) 
30 % given weights Wij, compute min dist Dij between node i to j 

% on shortest path from i to j, j has immeadiate predecessor Pij 

n = size(W,l); 
if(n~=size(W,2)) 
35 error('Input W is not square??! !'); 

end 

D = W; 

40 P = repmat([l:n]',[l n]); 

P = P .* ~isinf(W); 
P = P .* ~eye(n); 

for k = 1 ;n 
45 fori=l:n 
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end 

k; 

10 D; 

P; 
end 



for j = l:n 

alt_path = D(i,k) + D(k,j); 
if(D(io)>alt_path) 

D(ij) = alt_path; 

P(i,j) = P(k,j); 

end 

end 
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APPENDIX II 

function addlink(TOPO) 
% addlink(TOPO) 
% 

% interactively add links to the TOPO 

update(TOPO); 

c_src = 1 ; 
c_dst = 2; 
c_bw = 3; 

figure(TOPO.cur_fig) 
while (1) 

fprintf(i;\n\nHit Button 3 to end...\n\n'); 

% find coords and index i of src 
[xli yli button] = ginput(l); 
if (button == 3) break; end 

d = sqrt((TOPO.Iocs(:,l) - xli).^2 + (TOP0.1ocs(:,2) - yli).^2); 
[d,i] = min(d); 

xl =TOP0.1ocs(i,l);yl = TOP0.1ocs(i,2); 

% find coords and index j of dst 
[x2i y2i] = ginput(l); 

d = sqrt((TOP0.1ocs(:,l) - x2i).^2 + (TOP0.1ocs(:,2) - y2i).^2); 
[d,j] = min(d); 

x2 = TOP0.1ocs(j,l); y2 = TOPO.Iocs(j,2); 
hold on; 

Ih = line([xl x2],[yl y2],'color';red'); 
cap = input('Enter capacity (in Mbps) > '); 

:^rintf(l,' About to create symetric %d Mbps link fi:om node %d to node %d\n',cap,ij); 

doit = input('Enter Y to confirm, N to reject, and B to change bandwidth (Y)> ','s'); 

if (isempty(doit)) doit = 'Y'; end 

if(doit = 'n'|doit=='N') 
delete(lh); 
return; 
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end 

if (doit = 'b'|doit = ^B') 

buf = sprintfCEnter capacity from %d to %d (in Mbps) > ',ij); 
cap_i JoJ = input(buf); 

buf - sprintfCEnter capacity from %d to %d (in Mbps) > ' j,i); 
cap J_to_i = input(buf); 

else 

cap_i_to J = cap; 
capJ_to_i = cap; 

end 

%build the link records 
clear linkab linkba; 

linkab.src = i; 
linkab.dst = j; 
linkab.bw = capJ_toJ; 
linkab.handle = Ih; 

linkba. src = j; 
linkba.dst = i; 
linkba.bw = cap J_to_i; 
linkba.handle = Ih; 

% now draw the actual link on the map 
delete(lh); 

Ih = drawlink(TOPO, linkab); 

% now store the link info 

TOPO.links = [TOPO.links ; linkab ; linkba ]; 

TOPO.linkarray = [TOPO.linkarray ; [ i j capjjoj] ; [ j i capjjoj 

end % of while loop 
assigninCcaller',inputname( 1 ) JOPO); 



function [C, portmap] = capacity(TOPO) 
% [C, portmap] = capacity(TOPO) 
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% portmap maps indices of C to elts of nodes(TOPO) 
% [node dir] where 

% node is index of elt in nodes(TOPO) 

% dir is 1 if data enters here, -1 if data leaves here 

5 

numnodes = length(TOPO.Hnks) * 2; 

C = zeros(numnodes,numnodes); 

10 curnode = 0; 

portmap ~ []; 

for i = 1 :length(TOPO.Iinks) 
Iink = TOPOJinks(i); 
cumode ^ curnode + 1 ; 
15 portmap(cumode,:) = [link.src -1]; 

cumode = cumode + 1 ; 
portmap(cumode,:) ^ [link.dst 1]; 

C(curnode-l, cumode) = link.bw; 

20 end 

c_node = 1; 
c_dir = 2; 

25 for i = 1 :length(TOPO.nodes) 

ins = fmd(portmap(:,c_node) = i & portmap(:,c_dir) == 1); 
outs = fmd(portmap(;,c_node) == i & portmap(:,c_dir) == -1); 

for j = ins 

30 for k = outs 

CG,k)-inf; 

end 

end 

end 

35 function [a, b, c] = debug(t) 

update(t); 
fieldnames(t) 

40 

a = t.nodes 
b = t.locs 
c = t.links 

function display(TOPO) 
45 % DISPLAY a topo object 
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% a link is a unidirectional, so the value is probably twice what you want 

fprintfC[TOPO object: %d nodes %d Iinks]\n',... 

Iength(TOPO.nodes),length(TOPO Jinks)); 
5 function draw(TOPO) 

% draw(topo) 
% 

% draw the topology figure in a new window 

1 0 TOPO.cur_fig = figure; 

axis(TOPO.axis); 

axis equal; 

axis manual; 

box on; 
1 5 hold on; 

for i - l:length(TOPO.nodes) 

nm - plot(TOPO.nodes{i}.loc(l),TOPO.nodes{i}ioc(2);ob'); 
TOPO.nodes{i},mark_handle = nm; 
20 if(isfield(TOPO.nodes{i};nameloc')) 

TOPO.nodes{i}.nameloc(3) = text(TOPO,nodes{i}.nameIoc(l)5.,, 
TOPO.nodes{i}.nameloc(2)5TOPO,nodes{i}.name); 

end 

end 
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% yes, this draws the same link twice, fix it if it matters -dam 11/21 



T0P0.1inkarray=[]; 
for i = l:length(T0P0.1inks) 
30 T0P0.1inks(i).handle = drawlink(TOPO,TOPO.Hnks(i)); 

TOPO.linkarray = [TOPO.linkarray ; ... 

[ T0P0.1inks(i).src TOPO.Iinks(i).dst T0P0.1inks(i).bw]]; 

end 

3 5 assignin('callef ,inputname( 1 ),TOPO); 

function ex(t) 

t.nodes 

function labelnames(TOPO) 
40 % function labelnames(TOPO) 

% make it easy to label the nodes 

for i = 1 :length(TOPO.nodes) 

fprintfCPlace label for node %d "%s"\n',i,char(TOPO.nodes{i}.name)); 

45 
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origcolor = get(TOPO.nodes{i}.mark_handle,'color'); 
set(TOPO.nodes{i}.mark_handle;color',[l 0 0]); 

if(isfield(TOPO.nodes{i};nameloc')) 

good_x ^ TOPO.nodes{i}.nameloc(l); 
good_j = TOPO.nodes{i}.nameloc(2); 

end 
th=[]; 
while (1) 

fprintfCButton 1 to (re)place text, Button 3 to accept\n'); 

[x,y, button] = ginput(l); 

if (3 button) break; end 

if (-isempty(th)) delete(th); end 

th = text(x,y,TOPO.nodes{i}.name); 

good_x = x; good_y = y; 

end 

TOPO,nodes{i}.nameloc = [good_x, good_y, th]; 
set(TOPO.nodes{i}.mark_handle,'color',origcolor); 

end 

assignin('caller',inputname(l),TOPO);function names(TOPO) 
% NAMES the list of names of the nodes in the topo 

Q)rintf(TN[ode\t\tName\n'); 
for i = l:size(TOPO.names,l) 

fprintfC%d\t\t%s\n',i,TOPO.names{i}); 

end 

function [node] = nodes(TOPO) 

% function [node] - nodes(TOPO) 

% returns a cell array describing nodes in the TOPO 

node = TOPO.nodes; 

function [TOPO] = topo(TOPO) 

% [TOPO] = topo(TOPO) 

%% if input TOPO is 'inif, create a new topology 

% 

% newtopo = topo('init'); 
% 

% else add new nodes to TOPO 
% 

% nodes is a array of structs, one per node 

% link is a array of structs, one per link 

% a link is a unidirectional item, so there are probably twice 

% as many links as you'd expect. 
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if (nargin< 1) 

error('topo(TOPO) or topo("init") - not enough args'); 

end 

if (ischar(TOPO) & TOPO 'inif) 
clear TOPO 

TOPO.nodes - []; 
TOPOJinks-[]; 

TOPO. capacity = []; % now computed as needed 
TOPO.locs = []; % internal cache 
TOPO.Iinkarray = []; % internal cache 

f = figure; 
axis([0 75 0 50 ]); 
TOPO. axis = axis; 
TOPO.cur„fig-f; 
axis equal 
axis manual 
box on 
hold 

else 

figure(TOPO.cur_fig); 

end 



nodecount ^ length(TOPO.nodes); 

while (1) 

clear nodeinfo; 

fprintf(i;\n\nHit Button 3 to stop\n\n'); 

[xybut] = ginput(l); 

if (but = 3) break; end 

x = floor(x); y = floor(y); 

nm ^ plot(x,y,'ob'); 

name ^ inputCEnter name > Vs'); 

nodeinfo. loc = [ xy]; 
nodeinfo.mark_handle === nm; 
nodeinfo.name = cellstr(name); 
nodecount = nodecount + 1 ; 
TOPO.nodes{nodecount} = nodeinfo; 

end 
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ifCtopo' class(TOPO)) 

TOPO - class(TOPO;topo'); 

end 

5 

if (nargout = 0) 

assignin('caller'4nputname(l),T0P0); 

end 

function Ih drawlink(TOPO, link) 
10 % assumes TOPO.linkarray is already valid, and draws the position of 

% link line based on the number of links already present in linkarray 

c_src = 1 ; 
c__dst = 2; 
15 c_bw=3; 

i ^ link.src; 
j link.dst; 

20 xl - TOPO.nodes{i}.loc(l) 

yl =TOPO,nodes{i}.loc(2); 
x2 = TOPO.nodes{j}.loc(l) 
y2 = TOPO.nodes{j},loc(2) 

25 

if (isempty(TOPO.linkarray)) 
num_links = 0; 

else 

num_links = sum(T0P0.1inkarray(:,c_src) = i & TOPOdinkarray(:,c_dst) == j); 

30 end 

pattern- [0 1 -12-23 -3]*. 3; 

if (abs(xl - x2) > abs(yl - y2)) 
35 delta__x = 0; 

delta_y = pattem(num_links +1); 

else 

delta_x - pattem(num_links +1); 
delta__y = 0; 

40 end 

Ih = line([xl x2] + delta_x, [yl y2] + delta_y, 'color', 'black'); 
function update(TOPO) 

45 
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clear TOPO.locs; 

for i - 1 :length(TOPO.nodes) 

TOP0.1ocs(i,:) = TOPO.nodes{i}.loc 

end 

clear TOPO.linkarray; 

for i - 1 :length(T0P0.1inks) 

TOPO.linkarray = [TOPO.linkarray ; ... 

[ T0P0.1inks(i).src T0P0.1inks(i),dst T0P0.1inks(i).bw]]; 

end 

% these are here to be cut and pasted into other functions as needed 

% there doesn't seem to be a good way to pass them around in another fashion 

% (using assigning('caller\..) to force their definition sounds like asking 

% for trouble 'cause you'll overwrite another definition of them...) 

c_src - 1 ; 

c_dst = 2; 

c_bw = 3; 

assignin('caller',inputname(l),TOPO); 



